Density profiles for atomic quantum Hall states 
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We analyze density profiles for atomic quantum Hall states, which are expected to form in systems 
of rotating cold atoms in the high-rotation limit. For a two-dimensional (single- layer) system we 
predict a density landscape showing plateaus at quantized densities, signaling the formation of 
incompressible groundstates. For a set-up with parallel, coupled layers, we predict (i) at intermediate 
values of the inter-layer tunneling: a continuously varying density profile p(z) across the layers, 
showing cusps at specific positions, (ii) at small values for the tunneling, quantum Hall-Mott phases, 
with individual layers at sharply quantized particle number, and plateaus in the density profile p(z). 
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Among the fascinating developments in the field of 
quantum gases is the possibility to study correlated states 
of matter in a setting that is entirely different from the 
traditional setting of electrons in a solid state environ- 
ment. A prime example are fractional quantum Hall 
(qH) states, which are expected when trapped atoms 
(bosons or fermions) are made to rotate at ultra-high an- 
gular momentumi*2i2i2i£i£i£. While there has been steady 
progress in achieving high rotation rates 8 -, the conditions 
for the actual realization of these states have not yet been 
met. [Se o 9 i 10 for alternative proposals aimed at condi- 
tions where atomic qH states are expected.] 

The most direct experimental signature of electronic 
(fractional) qH states, the quantization of the Hall con- 
ductance, is not easily available for realizations of such 
states with neutral atoms. It is thus important to investi- 
gate the observable features of atomic qH states. Earlier 
studies have focused on the fractional statistics^, loss of 
condensate fractioni^, properties of edge excitations^, 
and on density correlations in expansion images^. 

In this Letter, we work out the proposals put forward 
iuiSjiJi f or t ne detection of atomic (fractional) qH states 
via characteristic density profiles. We present results for 
a single layer, and for multi-layer configurations. We 
shall focus on the fractional qH states formed of interact- 
ing bosons, although qualitatively similar results would 
be obtained for the fractional qH states of interacting 
fermions. We remark that non- interacting fermions in a 
rapid rotation regime, and subject to a slowly varying 
potential, will display density profiles similar to the ones 
we discuss herein. In that case, the effects arise from the 
Landau level structure for non-interacting fermions. All 
results discussed in the present work are strong many- 
body interaction effects. 

We first analyze a single layer situation, with rota- 
tion at or near the critical frequency lo±. The char- 



acteristic length scales £± and £u are set by the har- 
monic confinement lo±_ in the x — y (in-plane) direction 
and lou in the z (out-of-plane) direction, according to 
£x,\\ = \Jh/ ( mLU ±.\\)- The energy scale characterizing 
the qH states is g q n = ( 2 J)3/2 mi)|i° s w ^ n a s the scatter- 
ing length. Assuming a s — 5 nm, fkoj_ ~ 5 nK, £m = 50 
nm, (? q H is in the order of 1 nK. 

The state of matter formed by rapidly rotating bosonic 
atoms, at critical rotation lo — lo±_ and in the absence 
of any additional potential, depends on the filling frac- 
tion v = n/no [no = l/irijj^. For v less than v c ~ 6, 
the vortex lattice is unstabl e) 3 ! 18 , and the groundstates 
are homogeneous quantum fluids. Their interaction en- 
ergy density e[u] is a complicated function of v, contain- 
ing cusps that indicate incompressible groundstates at 
specific filling fractions v^. The incompressible ground- 
states have been studied by exact diagonalization studies 
in edge-less geometries (sphere or torus), where it was 
found that they are well-described by a variety of (frac- 
tional) qH liquidsiiSi 3 ^*^. These include the Laughlin 
state at v = 1/2, a series of states at = p/(j> + 1) 
that can be understood in terms of composite fermions 
(CF) 2 , and non-abelian states such as the Moore-Read 
(MR) and Read-Rezayi (RR) states^ at v = k/2. 

Here we focus on inhomogeneous qH liquids that form 
in the presence of a residual potential V(r) in the rotating 
frame of reference. We consider (i) rotation lo slightly be- 
low the critical frequency lo± , leaving a residual parabolic 
potential V^r) = ^k^r 1 with ki oc (lo± — lo), (ii) critical 
rotation lo = lo±, with an additional confining potential, 
here taken to be quartic, V±{r) — jk^r 4 . We shall dis- 
cuss first the case in which there is a very large number 
of atoms, N 3> 1, in which case mean density profiles can 
be reliably found from a one-shot measurement of parti- 
cle positions; we shall then turn to discuss the effects of 
fluctuations when N is not very large. 
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If the potential V(r) varies slowly in space compared 
to £±, the density distribution, n(r), averaged on scales 



large compared to £±, can be obtained by minimizing 

E = Jd 2 r{e[v(r)]+V{r)n(r)~iin{r)} , (1) 

where e[v] is the interaction energy density arising from 
the contact interactions, v{r) = n(r)/no and [i is the 
chemical potential. That the energy functional is lo- 
cal significantly simplifies its minimization: at each po- 
sition, the density n(r) is that which minimizes e[^(r)] — 
HLn(r), where the local chemical potential is [i^ (r) = 
fi — V(r). For a vortex lattice (at large v), the energy 
density is e\v\ = noi/ 2 b g q n, where b = 1.1596 is the 
Abrikosov parameter for a triangular lattice. The depen- 
dence of density on /ii(r) is then simply n(r) oc hl(t), 
leading to a Thomas-Fermi profiled For the qH regime, 
v < v c , there are cusps in the energy function e[v], which 
give rise to a step-like dependence of n on \xl ; this causes 
plateaus in the spatial density distribution. 

As a simple example, we consider a situation where 
liquids at v = 1/2 and v = 2/3 make up the lowest 
energy configuration. The energy per unit area takes the 
values e[v = 1/2] = 0, and e[v — 2/3] = (2/3)noe 2 (with 
e 2 a number of order g q n, see-). We find v = for fi < 0, 
v = 1/2 for < fi < 4e 2 , and v = 2/3 for fi > 4e 2 . 

Let us take a harmonic confining potential V(r) = 
^fc 2 r 2 , such that /Uc(r) = M — \kir 2 . The v = 1/2 state 
forms a disc that extends out to fii, = 0, i.e. to a ra- 
dius ri — y/2fi/k2. There is a critical value of chemical 
potential, \x c — 4e 2 , at which the v — 2/3 state will first 
appear in the center of the trap (where /i_L is maximum) . 
At p, = fj, c , the number of atoms in the v — 1/2 disc is 
N c = l/2n 7rr 2 = 2e 2 /A 2 , with A 2 = \k 2 i\. Above this 
critical value, the inner disc at v — 2/3 has a radius r 2 = 
y/2(fi — Hc)/k2- Expressed in terms of N and N c , the 
locations of the two steps are r 2 /£± — \/3(N — N c )/2, 
rx/lx. = y/(3N + N c )/2. From the dependence e[v] ob- 
tained from exact diagonalizations 3 it can be inferred 
that intermediate densities other than v = 0,1/2,2/3 
will not appear in the density profile. 

Repeating this analysis in a quartic potential T4(r) = 
r 4 leads to a critical N of N c = yia/A^, with 
A4 = jk^ 2 ^, with inner and outer edges at (r 2 /^) 2 = 

3 (v / 9A r2_ 87V 2-A0, ( n /£ ± f = iN-ym^-wi 

Including in the analysis more quantum liquids, at fill- 
ing V\ < i>2 < ■■■! leads to a sequence of critical val- 
ues N — Nc k \ marking the onset of the formation of a 
central region of a qH fluid at v = v^. To determine 
the corresponding density profiles we use the interaction 
energy function e[v\ calculated by exact diagonalization 
studies on a torus 3 . In the graphs presented in Fig. ^ 
a confining potential V(r) — |fc 2 r 2 is assumed. The 
energy functions used are for a system size of Ny = 6 
single-particle states on a torus with aspect ratio V3/2 
(solid lines) [with numerical results up to v = j], and 
for Ny — 12 single-particle states on a torus with aspect 
ratio 0.3 (dashed lines) [up to v — |]. For filling fractions 
larger than v c ~ 6, the groundstate is a compressible vor- 
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FIG. 1: Radial density profiles in a single layer with harmonic 
confinement, for a series of filling fractions at the center of the 
trap. The solid and dashed lines correspond to energy func- 
tions e[v] obtained from numerics using Ny = 6 and Nv ~ 12 
states, respectively. 



tex lattice^. In this regime, the density profile (averaged 
on a lengthscale large compared to £±) is an inverted 
parabola 20 . The steps seen in Fig. 1 in this regime are 
an artifact of the finite-size numerics (on a finite system 
with Ny = 6 the minimum change of filling fraction is 
Av = 1/6). For filling fractions less than v c ~ 6, plateaus 
appear at the filling fractions of incompressible qH fluids, 
including the sequence of MR and RR states at v = k/2 
[horizontal dashed lines] as well as some of the CF se- 
quence v — p/ (p+1) withp = 2, 3, —4 [dotted lines]. This 
plateau structure is a direct consequence of the incom- 
pressibility of the fractional quantum Hall liquids. This 
structure could be observed in the density distribution of 
a single-layer system following expansion 14 . 

We have tested the above picture against exact diag- 
onalization results for small systems. However, at the 
small sizes that can be handled in exact diagonalization 
(up to N — 12 particles on a disc), the density profiles 
seem to be dominated by edge effects, indicating that the 
system sizes do not exceed the correlation lengths in the 
expected qH liquids^ 1 -. 

We next investigate the situation where an optical lat- 
tice in the z direction is imposed on a cloud of atoms ro- 
tating around the z-axis^If the optical potential is suffi- 
ciently deep, this will define a stack of parallel planes with 
(weak) tunneling t between the planes. The lengthscale 
£11 is now determined by the confinement frequency in a 
single layer of the optical lattice. We assume that there is 
an overall quadratic confinement in the z-direction, giv- 
ing rise to a chemical potential fi(z) = fiQ — /i 2 (z/d) 2 
(with d the distance between the layers), which will in- 
duce a slow modulation in the number of particles per 
layer. The idea is that parameters can be chosen such 
that, while the original cloud is in a mean field regime 
(y > v c ), displaying a vortex lattice, the individual lay- 
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ers denned by the optical lattice can be in a quantum 
regime, so that the entire configuration becomes a stack 
of weakly coupled quantum liquids. With Nl layers, a 
total filling fraction v (for the entire cloud) gives a filling 
fraction per layer of v' = v/Nl- If v' < v c , the vortex lat- 
tice will melt if the inter-layer tunneling energy t is made 
smaller than a critical value t Cl . We have generalized the 
analysis oii& to determine the value of t Cl (seeSi for a 
detailed account). For N = 5000 particles, N v = 100 
vortices and a number of layers Nl — 50, our estimate is 

t ci ~ O.lffqH. 

Assuming N/Nl ^> 1 and t -C g q n we have the fol- 
lowing physical picture. Within each layer, we have a 
density landscape built out of incompressible qH liquids 
each of the form shown in Fig. 1. This landscape varies 
slowly from layer to layer. The gap for bulk excitations 
over each qH liquid is of order gqtr^ and since t <C <7 q n 
there will be no 'bulk-to-bulk' tunneling events between 
the layers. However, the energy scale for processes where 
atoms tunnel from the edge of a qH fluid in one layer i to 
the corresponding edge in an adjacent layer i± 1, is much 
lower, of order <7Mott- This scale is defined as <?Mott — A2 
for quadratic confinement and <7Mott — XiN for the quar- 
tic case, with N{ the number of particles in the qH liquid 
in layer i. In terms of N = 5qH/A 2 (N = v/5qH/A 4 ) 
we have g Mo tt ^ g q H/N (gMott ^ g q HNi/N(j) for the 
quadratic (quartic) case, if the tunneling t is large on the 
scale of gMottj inter-layer tunneling events will establish 
a continuously varying density profile p{z) — (iVj) across 
the layers (where N is now the total number of parti- 
cles in layer i). This profile could be measured from an 
image of the confined condensate, provided the period of 
the lattice, d, is made larger than the imaging resolution 
(perhaps by forming the optical potential from running 
waves at a shallow angle) M- 

We first assume quadratic in- plane confinement. If the 
number of particles per layer stays below Nc = Itij A2, 
the qH liquids in the layers will be Laughlin liquids at 
v = 1/2. Minimizing the total energy w.r.t. each of the 
Ni, we obtain a parabolic TF density profile 

with z\ — d 3 ^j-^. If N, H2 are such that some of the N 

(k) 

exceed the critical values Nc , there will be additional 
stable qH liquids in these layers. The corresponding pro- 
file p{z) exhibits cusps at the positions ±Zk where N goes 

Ik) 

through Nc ■ If the central layers near z = support 
v = 1/2 and v = 2/3 liquids, we have 

P ^ = TP2^^~ z ^ for Z2< \ z \ <Zl 

= ^ (2) + i^^ 2 - z2 ) fOT W<* (3) 

with z\-z\ = 2&N<M and z\ + z|/3 = d 3 ^-^. This 
profile shows a (weak) cusp at \z\ — Z2, with the slope 
■§-p{z) changing by a factor of 4/3. 
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FIG. 2: Density profiles p(z) for quadratic (top) and quartic 
(bottom) in-plane confinement. The layers at —Zi < z < —Z2 
have the v\ = 1/2 liquid only; the cusps at z — —zi mark the 
onset of a second liquid with V2 — 2/3 in the layers. 
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FIG. 3: Schematic density profile p(z), for quartic in-plane 
confinement and qH liquids at ui = 1/2, V2 = 2/3, vz = 1. 

The nature of the cusps at z — Zk depends on the 
in-plane potential V(r). In Fig. [3 we show the cusp at 
z = —Z2 for two choices of V(r). Clearly a steeper con- 
finement in the layers leads to a more pronounced cusp 
in the profile p(z). [A square well cut-off in V(r) can 

produce a flat plateau in p(r).] For quartic confinement, 

(2) 

the profile at N < Nc is a semi-circle, and the cusps at 
\z\ = Zk have square-root singularities. In fact, p(z) is of 

the general form p(z) = J2k a fcV ' z 1 z2 - Fig-El s h ows 
p(z) for quartic confinement, with the layers supporting 
quantum liquids at V\ = 1/2, v% = 2/3, ^3 = 1. (We note 
in passing that similar considerations lead to the conclu- 
sion that the density distribution following an expansion 
of a multi-layered system in an otherwise quadratically- 
confined trap will have a dependence on the radial co- 
ordinate r that is of the same qualitative form as Fig. 
is of the axial co-ordinate z.) 

We now turn to consider the situation when N is not 
very much larger than 1. We begin by estimating the 
fluctuations. A deviation 8N from the mean field profile 
costs an energy of order gMott{SN) 2 , giving 

SN^.In^, SN^J^-L (4) 
V 9qH V N g qH 
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FIG. 4: Schematic density-profile p(z), for layers supporting 
Laughlin liquids in quartic confinement. The horizontal seg- 
ments indicate the qH-Mott phases, with sharply quantized 
particle number per layer, while the remaining parts corre- 
spond to 'conducting' phases, with fluctuating Ni. 



for the quadratic and quartic cases. In addition there will 
be intra-layer fluctuations, for example when a particle at 
radius r$ is promoted to an unoccupied orbital at fj +5r, : . 
Equating the corresponding energy to t leads to 



(Si 



(Sn 



No 
N 



t 

.9qH 



(5) 



for the quadratic and quartic cases. This shows that if 
t *C 5 q H and Ni/No is not too small, the steps in the 
in-layer density profiles are well-defined on the scale of 
the overall radius. 

A one shot experiment will obviously produce integer 
values for all N{. Based on the above, we predict that 
these numbers will follow our continuous curves for p(z) 



rather closely, with fluctuations as specified in Eq. J3J. 
Gradually lowering t (for example, by turning on the po- 
tential slowly enough) will lead to values Ni that are 
nearest integers to the corresponding values of p(z). 

The estimates in Eq. (0J indicate that 6Ni ~ 1 if t is 
lowered to a value of order gMott- Based on the analogy 
with single atoms in an optical lattice, we anticipate a 
Mott phase where the numbers Ni are sharply quantized. 
We have analyzed an effective Bose-Hubbard model, with 
layers supporting a v = 1/2 state in quadratic confine- 
ment. In a mean-field treatment we find, for fj, satisfying 
2A2(iV — 1) < fi < 2X2N, the following critical t for en- 
tering the TV-particle qH-Mott phase 



t C2 (N) = -\ 2 [N ■ 



2A 



-][{N-1) 



JL] 

2A 2 J 



(6) 



This has a maximum of A2/4, which is of order ^Mottj 
for all N. For confinement 'steeper than' quadratic, the 
maximal t C2 (N) grows with N, consistent with gMott — 
X4N for the quartic case. 

Scanning the sample from layer to layer, there will then 
be segments showing qH-Mott phases, and 'conducting' 
regions in between, see Fig.0] A measurement of the par- 
ticle number per layer in an A^-particle qH-Mott phase 
will give Ni = N with no fluctuations. The same mea- 
surement in a 'conducting' phase will produce integer 
numbers as well, but these will show fluctuations. 
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